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Abstract. We prove global and local upper bounds for the Hessian of log positive 
solutions of the heat equation on a Riemannian manifold. The metric is either fixed or 
evolves under the Ricci flow. These upper bounds supplement the well-known global 
lower bound. 



1. Introduction 

Gradient bounds for log solutions of the heat equation have appeared in the important 
papers by Li and Yau [II] and by Hamilton [7J. The main result in |llj can be regarded 
as a lower bound of the Laplacian of log solutions under the assumption that the Ricci 
curvature is bounded from below. These bounds came in both global and local versions. 
The main result in [7] is a global lower bound for the Hessian matrices of log solutions 
under certain curvature assumptions. Many applications of these results have been found 
in numerous situations. See for example the papers [I], [2], [3], [I], [5], [5], [ID], [H], [13] . 
[H] and also the books [9J and [15] . 

In this paper, we derive an upper bound of the Hessian matrices of log solutions of 
the heat equation on manifolds. We prove both a global and local version of the bounds 
which take two different forms and which are generally sharp in respective cases. While 
the global version can be proven by building on the ideas in [7] , the local version requires 
additional quantities and calculations and appears to be unknown even in the Euclidean 
setting. We give a unified proof for both versions. In addition, we generalize the Hessian 
bound for log solutions to the Ricci flow case. Interestingly, the Ricci flow induces a 
cancelation effect which makes the curvature assumption less restrictive than the fixed 
metric case. 

In order to present results, let us fix some notations to be used throughout the paper. 
We denote by M a Riemannian manifold with metric g. For simplicity of presentation, 
we always assume that M is a compact manifold without boundary in this paper. The 
local bounds clearly hold without this assumption. The global bounds also hold when 
one imposes suitable conditions at infinity. We use Ric and Rm to denote the Ricci and 
full curvature tensor respectively. In local coordinates, the metric g is denoted by gij, the 
Ricci curvature by Rij etc and the curvature tensor by Rijki etc. We use the convention 
that R^ = g kl Ritij in local coordinates. The Hessian of a function u is written as Uij. 
If V is a 2-form on M and £ is a vector field, then in local coordinates, we use £ T V£ to 
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denote V (£,£). The distance function is denoted by d(x,y) in the fixed metric case and 
by d(x,y,t) in the Ricci flow case. A geodesic ball is denoted by B(x,r) or B(x,r,t) 
where x is a point in M and r is the radius, and t is the time when the metric changes. 
For R > and T > 0, a parabolic cube is defined by 

QR,T(x ,t ) = {(x,t) | d(x ,x,t) < R,t -T <t< t }. 

For the fixed metric, we simply have 

Qr,t(xq, t ) = B(x , R) x (t - T, t \. 

A positive constant is denoted by c and C with or without index, which may change 
from line to line. 

We first state our upper bound of the Hessian matrix of positive solutions for a fixed 
metric. 

Theorem 1.1. Let M be a Riemannian manifold with a metric g. 
(a) Suppose u is a solution of 

d t u - An = in M x (0,T]. 

Assume < u < A. Then, 



t(uij) < u(5 + Bt) 1^1 + log —J inM x (0,T], 

where B is a nonnegative constant depending only on the L°° -norms of curvature tensors 
and the gradient of Ricci curvatures. 
(b) Suppose u is a solution of 

d t u-Au = in Q RtT (x ,to)- 

Assume < u < A. Then, 

(uij) < Cu + + b\ fl + log^ in QR T(x ,to). 

where C is a universal constant and B is a nonnegative constant depending only on the 
L x '-norms of curvature tensors and the gradient of Ricci curvatures in Qr^t{xo, to). 

In Theorem ll.il Part (a) is the global version and Part (b) the local version. We point 
out that different powers of l+log(^4/w) appear in the right-hand sides. We will illustrate 
by an example that the extra power in the local version is optimal. This phenomenon 
adds more variety to the long list of differential Harnack inequalities. 

We remark that the global estimate in Part (a) can be obtained as a consequence of 
the lower Hessian estimate and the upper Laplace estimate by Hamilton [7j. Arguments 
in this paper work for both global and local estimates, and they extend to the Ricci flow 
case. 

Next, we state our upper bound for the heat equation coupled with the Ricci flow. It 
has a similar form to the fixed metric case. 
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Theorem 1.2. Let M be a Riemannian manifold with a family of metrics g = g(t) 
satisfying 

d t g = -2Ric in M x (0,T). 

(a) Suppose u is a solution of 

d t u - An = in M x (0,T]. 

Assume < u < A. Then, 

t( Uij ) < u(18 + Bt) ^1 + log ~\ in M x (0, T], 

where B is a nonnegative constant depending only on the L°° -norms of curvature tensors. 

(b) Suppose u is a solution of 

d t u - Au = in Qr,t(xq, t ). 

Assume < u < A. Then, 

(uij) < Cu + + b\ ^1 + log^ in QR T(x ,t ). 

where C is a universal constant and B is a nonnegative constant depending only on the 
L 00 -norms of curvature tensors in QR,T(%o,to). 

We point out that the constant B in Theorem 11.21 does not depend on the gradient of 
Ricci curvatures. 

We now describe briefly the method to prove both results. As expected, the general 
idea of the proof is still the Bernstein technique of finding a quantity (auxiliary function) 
involving the Hessian, which satisfies a nonlinear differential equation amenable to the 
maximum principle. The work is to find such a quantity. In this paper, the central 
quantity is the quotient of the Hessian of the solution with the density of Boltzmannn 
entropy of the solution, i.e., , where u is a positive solution of the heat equation. 

Another quantity is jnn« - A Bernstein type argument is carried out for a combination 
of these two quantities and suitable cutoff functions as in [TT] . 

The rest of the paper is organized as follows. In Section [21 we prove Theorem 11.11 
about the Hessian bound on a Riemannian manifold with a fixed metric. In Section [31 
we prove Theorem 11.21 which treats the heat equation coupled with the Ricci flow. 

Acknowledgment. We wish to thank Professor Lei Ni for helpful suggestions. 

2. Heat Equations under the Fixed Metric 

Let M be a Riemannian manifold with metric g and A be the Laplace-Beltrami 
operator. We consider a positive solution u of the heat equation 

u t = Au in M x (0, oo). 

We assume 

< u < A. 
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Set 

f 1 U 

/ = log T 

Let {x 1 , ...,x n } be a local orthonormal frame at a point, say p 6 M. Then 

/. „ ity UiUj 

u u u z 

and hence 

ft = A/ + |V/| 2 . 
We first derive two equations on which the theorems are based. 

Lemma 2.1. Set 

Uij 

Then, 

-d t + A - JLvf ■ v) % - = 

+ ^ - {-2R kij iUki + RuUji + iijjiitf + (Vi-Rj7 + VjRu - ViRij)ui\. 
Proof. By noting 

dt{u(l - /)) = -u t f, d k (u(l - /)) = -u k f, 

we have 

dv— — + Ui ^ n * 



„(!_/) U 2(l-/)2' 
OkVij = —r 7T + 



u(l-/) u2(l-/)2- 

Recall the commutation formula (see [5] p219 e.g.): if Am — dtu = 0, then the Hessian 
satisfies 

-d t Uij + Auij = -2R kij iu k i + RuUji + + (ViRji + VjR a - V/i?^)^. 

Note 



_ Au^ UjjfAu 2u ijk u k f Ujjf k u k 2uijf 2 u 2 k 

' Hj ~ U(l -f) + U 2 (l - /) 2 + _ /)2 + U 2 (1 _ /)2 + U 3 (1 _ /)3 ' 



With tifc = u/fc and the commutation formula, we then have 
(A - d t ) Vij = ^4 + J^jWP + **f\ Vf \* 



+ u ^ — [-2R k ijiUki + RiWji + Rjiuu + (ViRji + Vjiij; - X7iRij)ui] 



2ffi 



k 



1-/ 



^ijfc Uijff k 



u(l-f) u(l-f) 



, ^ iv/l 2 

(1-/) 



+ ^tj- — -p[—2R k ijiu k i + RuUji + + {ViRji + VjRa — ViRij)ur]. 
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With the help of the expression for d k Vij, we obtain 
(A - d t )Vij = - -JkOkVij + Vij 



+ u ^ _ m [-ZRkijiUki + RuUji + Rjiuu + (ViRji + VjRu - ViRij)ui\. 

This is the desired result. □ 
Lemma 2.2. Set 

UiUj 



10 « 2 (1-/) 2- 

T/ien, 

-9t + A — -^7 V / • V ) 
2|V/| 2 

= — ^ + 2(t, w + f Wkl )(v kj + f Whj ) + R tkWk] + R 3M 

Proof. We proceed similarly as in the proof of Lemma l2.1i First, 

uuUj+UiUtj 2u t UiUjf 
(htfjj = — ^f- J-f- H 



u 2 (l-/) 2 u 3 ( i_ /)3 ' 

UfcjUj + UjUjfc 2u k UiUjf 
dkWij = — 7a2 + 



n 2 (l-/) 2 «3(l-/)3- 

Using Bochner's formula, we arrive at, after differentiation, 

(AujiUj + 2u ki u k j +Ui(Au)j u k uj u k Ui 

*J ~~ 2?1 ?A2 27l ?A2 ' 3 k ~ 



V?{l-ff l %2 (1 _ /)2 ^ U 2 (1 _ /)2 



^{u ki Uj + u kj Uj)u k f 2ujUjAuf 2ujUjU k f k §UiUjU 2 k f 2 



Hence 
where 



U 3(l-/)3 n3(l-/)3 n3(l-/)3 u\l - f)^ 

(A - 9f)tyy = H + -Rjfc^fcj + RjkWki, 
_ 2u k jUkj _ 4:(u k iUj + UkjUi)u k f 2ujUjU k fk 6uiUjU 2 k f 2 

* "* O/-. /*\0 ~T~ Q/-. /»\ Q ~T~ Q/-i /'XQ "T 



U 2 (l-/) 2 u3(l-/)3 ^(l-/)3 u4(l-/)4- 

In the expression of iJ, we write 4 = 2 + 2,6 = 4 + 2 and hence 
_ 2(u ki Uj + UkjUi)u k f AujUjulf 2 2ujUjU k f k 

U 3(l-/)3 + u 4 (1 _ /) 4 + u 3 (1 _ /) 3 

2u ki u kj 2(u ki Uj + u kj Uj)u k f 2ujUjulf 2 

+ U 2 (l-f) 2+ u3(l-/)3 

2ttfc/ fukiUj+UiUjk 2u k UiUjf \ 2u k fk UiUj 



u(l-f)\ u 2 {l-f) 2 v?{l-ff) n(l-/) u 2 (l-/) 

, U ki + UjUfc/ \ / U kj ^ UjUkf 



u(l-f) u 2 {l-f) 2 ) WW) n 2 (l-/) 
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With Uk = ufk, the expression of dkWij and definitions of Vij and u>ij, we have 

H = u (l-f) dkWij + ^-f Wij + 2(yVki + f Wki )( Vk i + /^fcJ')- 
This implies the desired result. □ 
Remark 2.3. We define the trace w of (w{j) by 

iVnl 2 



w = tr(wij) 



« 2 (1-/) 2 ' 
We also have 

-d t + A - y^y V/ • V ) %' = C 1 " f) wv v 

+ u ^ - [-2R kij iu k i + RiiUji + Rjiuu + (ViRji + VjiZj/ - ViRi^ui], 

-dt + A- -^rjVf ■ = 2(1 - f)w Wij 

+ 2(^fci + fw k i){v k j + fw k j) + RikWkj + RjkWki- 
We now prove Theorem 11.11 



Proof of Theorem \l.l\ Part (a). We first perform some important calculations. Let p 
be a point on the manifold and {x 1 , x n } be a local orthonormal coordinates system. In 
this coordinate and using the same notations as in Lemma [2. II and Lemma [2.21 the (2, ex- 
tensor fields and iOy can be regarded asnxn matrices. Set V = (%), W = (w^), 
w = tr(W), and 

L = -d t + A - r^-jVf • V. 
Then, by Lemma 12.11 and Lemma 12.21 

(2.1) LV = (1 - f)wV + P, 

(2.2) LW = 2(1 - f)wW + 2(V + fW) 2 + Q, 
where P and Q are matrices whose (i, j)-th components are 

P*i = ~~[\ 7\[—2RkijiUki + RuUji + Rjiuu + {V{Rji + VjRu — ViRij)ur] 

= —2RkijiVki + Ru v ji + Rjivu + (yi.Rji + VjRu — Vfit!ij) ^ — , 

and 

(2.4) = RikWkj + RjkWki- 

For a constant a 6 R to be determined, we have 

L(aV + W) = a(l - f)wV + 2(1 - f)wW + 2(V + fW) 2 + aP + Q. 
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Let £ 6 TpM be a unit tangent vector at the point p. We use parallel translation along 
geodesies emanating from p to extend £ to a smooth vector field in the local coordinate 
neighborhood. We still denote the vector field by £. Since V and W are (2, 0)-tensor 
fields, the function 

A = ^(aV + W)£ = {aV + W)(£, 
is a well-defined smooth function in a neighborhood of p. Then 

LX = H + f(aP + Q)£, 

where 

(2.5) = a(l - /)< T F£ + 2(1 - /)< T W£ + 2|(y + jW)£| 2 . 

By a£ T V£ = A - i' ll i. we have 

H = (l- f)w(\ - ewt) + 2(1 - f)w?W£ + 2\{V + fW)i\ 2 

= (i - f)w\ + (i - /k t w£ + 2\(v + /^)ei 2 . 

To simplify the last term further, we fix the point p and assume £ is the vector field 
generated, via parallel translation through geodesies emanating from p, by an eigenvector 
of aV + W at p, i.e., at p, 

(aV + W)£ = A£. 

Then 

(y + f W )i = -£ - - W£ + /W£ = -£ - (- - f) W£, 
a a a \ a J 



and hence 



Hence 



\(v + fww = ^--[z-f)ewz+(±-A \wz\ 2 . 



A 2 






a 2 




V a 



# = + A ( w ~ ^ Tw t) ~ f X ( w ~ -? W t 



+ {l-f) w fWZ + 2 --/ 



a 



The last two terms are independent of A and nonnegative. Hence 



(2.6) H>^- + \(w- -U T M^) -f\(w- -( r Wi 



For the last term, we note W is a rank-one matrix and hence 

< w. 

With / < and choosing a > 4, the third term in nonnegative, if A > 0. If A > 0, the 
second term is also nonnegative. Hence 
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In summary, if A > 0, then at the point p, 

o\2 

(2.8) LX > —y- + £ T (aP + 

or 

Now we proceed to prove Part (a). Let t be a universal constant to be fixed later. 
With a = 4, suppose the 2 form 

aV + W - -g 

assumes its largest nonnegative eigenvalue at a space-time point (pi,£i), with t% > 0. 
Let £ be a unit eigenvector at p\. We use parallel translation along geodesies emanating 
from p\ to extend £ to a smooth vector field which is still denoted by £. Set, in a local 
coordinate 

(2.9) /i = fiaV + W - Lg)£. 
and 

(2.10) A = f{aV + W)£. 

Then, both fi and A are smooth functions in a space time neighborhood of (xi,ti). Also 

t\ T , r TT r T 



Ln = L (A - -J = LA - -2 = H - -2 + £ J (aP + Q)£. 

Here P is given by (|2.5p . We now evaluate at (pi, ii). Since A — r/i has its nonnegative 
maximum at (pi,ii), we have, by 



> I ( A - 1) > ^ - 1 - \f(aP + Q)£| at ( Pl , t x ), 



or 



2A r 

(2.11) < ^ + |£ T (aP + Q)£| at(pi,ti). 

In order to bound [i and A from above, we need to find an upper bound for |£ T (ctP + (3)£| 
at (pi,ti). 

Let £ = (£i, ...,£„). By Q and we obtain 

|£ T (aP + Q)£|<K r P£| + |£ T Q£| 



< 



-2R kij iv k i + RuVji + Rjivu + (ViRji + VjRn - V/P; 



+ \£,i£,j{RikWkj + RjkWki)\ 
< |Ci^Q [-2P fcij/ ^ + P^ 7 + P^]| + C|VPzc|v^ + C|Pic||W|. 
Writing a^; = av k i + u>&/ — etc in the last line, we deduce 
|e T (aP + Q)£| 

(2.12) < [-2R kij i(av k i + w fci ) + Ru(avji + ii;^) + Rji(av u + w^)]| 



+ [-2fl fcii ,ti; J H + flatly + Pj7^]| + C|VPic| vW + C|Pic||W|. 
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At the point pi, we can choose a coordinate system so that the matrix aV + W is 
diagonal. Let /xi, • • ■ , /x n be the eigenvalues of the matrix aV + W — jl, listed in the 
increasing order. Without loss of generality, we assume ui < and n n > 0. Then 

T T 

\Rkiji(avki + u>m)| < \Rkiji(avki + u>m - -4z)| + \Rkijihi\~ 

n 

< l^Rkijkiavkk + Wkk ~ -)| + C|i2m|- < C\Rm\(^ n + + C\Rm\-. 

k=l 



Similarly 

Combining the last few inequalities, we deduce 



T 

Ru{auVji + Wji)\ < C\Ric\(n n + |^i|) + C\Ric\-. 



\£ T (aP + Q)£| < C|i?m|(// n + | /^i |) + C\VRic\ % /\W\ + C\Rm\\W\ + C\Rm\j. 
In the following, we set K\ = \Rm\i,°o and K 2 = |V Ric\L°° • Then 

\t T (aP + Q)t\ < CK x {n n + \m\ + ^) + CK 2S /\W\ + CK x \W\ 
< CK x {ii n + l/ixl + 1) + C7X 2 + C(Kt + K 2 )\W\. 

Observe that 

m + (n- l)/x n > fj,i H h /in 

arUjj itjitj r \ aAu r 



Hence, 

aAu r 

"1 < (n - l)A*n 7Z 7T + n~. 

u(l-f) t 

Therefore 



|^ (aP + < C/Xi //„ - — + - + CiT 2 + C(ifi + K 2 )\W\ 

u(l-f) t 1 



By [TT], we have 

|Vd 2 Uf ci 

u z u t 

where Ric > —Kq for some Kq > 0. With = An, we get 

An ci 

< ^ + ^0- 

n i 

Since < n < A, we have 

1 



< 1. 



1-/ 

Therefore, at (pi,ii), it holds 

\F(aP + Q)£\ < C/Xx L n + iil ) -j- ( '( /v i Y, + /vY) + ( Vv) + Y.)|n'|. 
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By the definition of W = (wij), we have 



\W\< 



\Vu 



2 



By Theorem 1.1 in [7], we obtain 



|ir £ c ( - + k 



and hence 



|£ T (aP + Q)£| £ Cf^i + K 2 ) ^ n + i±Ij + C{K X K G + K 2 + K 2 K ). 

Since \i = [i n < A at (pi,ii), this shows, by (|2.1ip . that 

9\2 T / 1 i T \ 

< -z + C{K X + K 2 ) (A + + C^x^o + K 2 + K 2 K ) at { Pl ,h). 

A simple application of the Cauchy inequality yields 

A y/r + 1 „ 

- £ — + 5 at (pi,ti , 

a t 

where B is a nonnegative constant depending only on Kq,Ki and K 2 with the property 
that B = if K = K\ = K 2 = 0. Then 

T 1 

A - - £ (aVr + 1 -t)- +aB < aB at (pi,*i), 

by choosing r sufficiently large. In fact, for a = 4, we can take r = 8 + 4v5- 

By definition, /i = A — ? at (pi , fi) is the largest eigenvalue of the 2 form aV + W — jg 
on M x (0, T]. Therefore, given any unit tangent vector r/ S T X M, x G M, it holds 

/(^ + W)r? - ^(t?, /?) < (A - ~)|(p 1)tl) <«B in M x (0, T). 

Thus 

ir/^Vr? £ - + Bt. 
a 

This proves part (a) of the theorem. 

Part (b). Now we localize the result in part (a). It is unexpected that the local 
bound is different from the global one. We point it out in the remark below that the 
local form is also sharp in general. 

Let tp be a cutoff function which will be specified later. Then, for any smooth function 
rj, we have 

dt(i>r)) = d t i> r) + ipd t r], 
V(^ry) = T] + i/jS/t], 
A(iprj) = Aip rj + 2VV>Vr/ + ipArj. 
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Hence 

2/ 

ipLr) = —ip dtf] + ipAr/ — ip- — y^f ' 

2/ 

= -(d t {ipv) ~ Vdttp) + A(V>?7) - rjAip - 2Vip • V77 - :j- J -^V/(V(V>r?) - 77VVO 

2/ 

= -d t (iprj) + A(tprj) - - - V/ • V(V>t?) + ??<9 t ^ - ??A^ 

+ YZy^f ■ VV - 2VV • Vr/. 
For the last term, we write 

Vip ■ Vr] = —^-ipVrj = —^-iy (ijjrf) — Vtp rj) 
W> |V^| 2 

= — vw - — r?. 

Hence 

2f 2W 
V'Lr? = -d t {iprj) + A(ip V ) - ~^Vf ■ V(iprj) p- V(^) 

2/ 2|V^| 2 

+ r]dt4> - V^ip + rj- V/ • Vip H ; — r/. 

1 - / t/> 

Set 

(2.13) Ll = _a < + A --^-V/-V-^ V. 

Then 



^Lr? = Li(r]ip) - rjLxip, 



or 



Li(r]ip) = ipLr] + rjLiip. 
With A introduced before in (|2.10p . we have 

(2.14) LiOM) = ipLX + ALiV = + Z T {aP + Q)£] + XLiip. 

Here -ff is given by (|2.5p . Now we analyze L\ip. We write 

LiV = + A^ - ^£ - -^V/ • V^. 

The first three terms are obviously bounded by choosing suitable tp. For the last one, 
we write 

2/ „, 2/ r-^t Vl/J 



-/ 

Note < 1 and — = is bounded. We need to control 

1 - J vi> 
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To this end, we recall the equation for / 

-d t f + Af = -\Vf\ 2 . 

Then 

Lf = -d t f + A/ - Y^ylV/l 2 = -\Vff - ^- |V/| 2 = =l^l\Vf\ 2 . 

Note 

T^T - \ if/ ^" 3 - 

Then 

Lf>\\Vf\ 2 . 



Hence, we obtain 



T 2V^ 2W r~i\i f 



and then 



4 ri " V 
Now we consider, for some constant /3 G R + to be determined, 

iftLityX + /3/) = ^ 2 H + ^ 2 £ T (aP + Q)£ + V-ALiV + ^i/. 
In the following, we consider eigenvalues of the 2 form 

^(aV + W) + 

If £ is an eigenvector of ip(aV + W) + (3fg at some point (ar, i), then 

M<*V + W) + (3fg}£ = tf, 

or in local coordinates, 

^(aV + W)i = { l i-Pf)i. 
If if)(x, t) / 0, then £ is also an eigenvector of aV + W. Hence 

H = + pf. 

Again we extend £ to a vector field around x by parallel transporting along geod 
starting from x. The vector field is still denoted by £. 
Let be the parabolic cube given by 

n = Q R>T {x ,t ) = B{x ,R) x (t -T,t ]. 

Let 

I 1 = f[^{aV + W)+ Pfg]£ = + 0f. 

Then 
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We will estimate \i from above. Recall from (|2.14p and fi = A + fif that 
(2.16) ^L lt jL = i) 2 H + ip 2 Z T (aP + Q)£ + iPXL^ + ^L x f . 

We first have, from (|2. T[> . that 



cx 



o- 



V> z # > -^(^A) 2 + ^A ( w - -^i r Wi ) - /V^A ( w - -fWi 



a 



At points where > 0, we have 
and hence ip\ > 0. Then 



^A + /3/>0, 



a" 



By this, (I2TT51) and (f27T6|> . we deduce, 



Q- 



i W l + l^ + 35*!! + V?|v/i-^' 



V>A. 



For the last term, we use the Cauchy inequality to control the ^A factor by the first term 
to get 



- C \dti/>\ + |AV| + 



2|VV| 



2\ 2 



We now take 
(2.17) 

If c is sufficiently large, we have 



[3 = csup 



cy>|v/| 



2 IWI 2 



> -Csup 



v> 2 



and hence 



^i/i > A(V^A) 2 + V 2 e T («^P + Q)£ + ~/W/| 2 



c 



|^| + |A^| + 



^ IWI 4 

C sup — 



Let V' be a cut off function supported in the space-time cube Qr,t(xo, to) such that 
t/} = 1 in the cube of half the size QR/2,T/2( x Q,to). We also require that 



|VVI<|, |A^|<C^tl, 



14 
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Here Kq is a bound on \Ric\ as before. Also the quotients are regarded as when ip = 
somewhere. Without loss of generality we just take to = T. We also require that ip is 
supported in the slightly shorter space time cube Qr$t/±( x 0i *o)- As usual, the cutoff 
function can be constructed from the distance function, which is not always smooth. 
One can either mollify the distance by convoluting with a smooth kernel or use the well 
known trick by Calabi to get around singular points of the distance. 

Let hq be the maximal eigenvalues of ip(aV + W) + f3fg with a unit eigenvector £. 
Assume no 1S taken at the space-time point (pijii). Again, by parallel translation, we 
extend £ to a vector field in a neighborhood of p±, which is still denoted by £. We are 
interested only in the case hq > 0. Since / < and ip = on the parabolic boundary 
of fi, we know that p\ must lie in the interior of B(p, R). Define a function \x = /j,(x, t) 
around (pi, ti) by 

H = i T (iP(aV + W) + Pfg)Z. 
Since (pi,ii) is a maximum point of /x, we have 



> ^)L xl x > \{ip\f + V 2 e T («P + Q)£ 
or 



C 



|^| + |A^| + 



V? J 



C sup 



V> 2 



Hence, at (pi,ti), it holds 



(2.18) 



Prom (l2~T2l . 

^|£ T (aP + Q)e| 

< [-2R ki ji^{av k i + Wfci) + Rwip(avji + u^/) + Rjiip(avu + wu)]\ 

+ V [-2R kijl w kl + RnWji + RjiwuW + CiP\Ric\\W\ + Cty|V2Hc|v1W1 

< [—2Rkijiip(av k i + Wki) + R u ip(avji + u;^) + Rjiip(azva + u^)]| 
+ C(if + Ki + K 2 )iP\W\ + CK 2 . 

By splitting off a term f3f, we obtain 

V|£ T (aP + Q)£| 

< [ - 2R kij i{ip(av k i + Wfci) - f3f5 ki } + Ru{ip(avji + u^/) - /3/£ w } 

+ Rji{ip(av u + wu) - I3f5 kl }} 

+ C(K + Ki)^|/| + C(K + #i + #2,^1 + C^ 2 . 

Let /xi, ...,/x n be the eigenvalues of the 2 form ip(aV + W) + /3/<7 at (pi,ii), which are 
listed in increasing order. We assume without loss of generality that /xi < 0. Then the 
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above inequality implies 

4>\f(aP + Q)C\ < C(K + Ki)On + M) + C(K + K^^f] 
+ C(K + K x + if 2 )^|^| + Cty#2- 

Observe that 

fj,i + (n - l)/i n > /^i H h 



tr 



auij 



oeAu 
> 7T + 

Hence, 

. , aAn 

A*l < (« - l)Mn - V> — 7Z 7T + ™P / • 

«(! - /) 

Therefore, 



(2.19) ris v ~ — r «(l-/), 

+ C(Kx + K 2 )^|W| + CK 2 + CK x p\f\. 

Since Ric > — Kq, by [11] and our choice of the cutoff function ip, we have, for any 
a > 1, 

o/|Vd 2 uA /l 1 

We note that the time 1/T is actually 1/i in [11]. However, since our cutoff function 
is supported in a shorter cube, these two terms are equivalent. With ut = Au and 

(j, = tpX + fif < ipX, we get, at (pi,h), 

4> 2 \f{aP + Q)Z\ < CK^X + CK X (± + -L 

+ C(i^ + K 2 ) + CK 2 V> 2 |VF| + CKtPlfl 
From the definition of W = (iVij), we have 

|Vn| 2 
\W\ < — 

By Theorem 1.1 in |16j . we obtain 

4> 2 \W\<c(^ + -^ + K c 



and hence 



T R 2 



+ C{K X K Q + K 2 + K 2 K ) + CK x p\f\. 
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Substituting this to (I2.18p . we find, at (pi,£i), that 

1 ,(VA) 2 < CK^ 2 \ + C(tfi + K 2 ) f I + -L ) + r* ( - + 



a 2Vr ; " ir v 1 \T R 2 J \T R 2 

+ CiKtKo + K 2 + K 2 K ) + CKrflf], 

and hence 



where B is a nonnegative constant depending only on Kq,Ki and K 2 with the property 
that B = if K = K x = K 2 = 0. This implies 



^X + Pf<c(^ + ^ + B^+Cy/K^\f\+Pf. 
Since / < 0, we know that Cy/KM\ + 0f < CK X . Thus 

Mo = H^) = «A + /3/)| (pi , tl) < C + ^ + B) • 
Here £> may have changed from the last line. Therefore 

Hence, for any unit tangent vector £ at x with (x,t) S Qr^t, it holds 
^ T (a^ + ^ + /3/<C^ + -^ + ^ inQ^T, 

or 

j. / „ „, s . . „ / 1 1 



^ («^ + W)£ < C ( - + ^ + B ) + /3|/|, in Q R}T . 



Recall from (|2.17|) that (3 = j^p, we then have 



hT i I 1 



and hence 



4>e V£<C[~ + w +B)(l-f), 



II 



T R 2 



This implies the desired estimate. □ 

Remark 2.4. When we compare the local version with the global version, we note an 

A 

extra power of 1 + log — in Part (b) of Theorem 11.11 This turns out to be optimal. 

u 

Consider x = 2, R = 1, t = 2, T = 1 and Qi,i(2,2) = [1,3] x [1,2] C R x x R t . For 
a > 0, set 



u(rc, t) = e 



ax+a 2 t 
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This is a positive solution of the heat equation in Qn(2, 2). Note 

u xx{2,2) ^ 2 

and 



3a+2a 2 ■ ^ e 



3a+2a 2 



A= sup u = e da+ " a , log — — = log — - ^ = a. 

Qi,l(2,2) U{2,2) e 2a+2a 



TT ^:r:r(2,2) 

Hence, — ; — ana 

u(2,2) 



log- 



have the same order in a. 



u(2,2) 

3. Heat Equations under Ricci Flow 




In this section we consider the heat equation coupled with the Ricci flow on a manifold 
M, over a time interval (0, T], 

(3.1) 

The heat equation and its conjugate have served as a fundamental tool in the theory 
of Ricci flow developed by Hamilton and Perelman. The two authors and others have 
derived gradient estimate for positive solutions of the heat and conjugate heat equation. 
See the paper [5J, [J3], [TH], [3], [3] and [JJ for instance. Here we prove an upper bound 
on the Hessian of the log solution, which seems to be missing as the fixed metric case. 
The general idea of the proof is similar to that in the previous section. However, since 
the metric is evolving, there will be extra terms to deal with, especially the term RijUij. 
To treat this term, we need to use the latest result in pp which provides a Li-Yau type 
gradient bound in the Ricci flow case. 

We will keep the same notations as in the last section. First let us derive the evolution 
equation for (vij) and (iVij). In this situation, the corresponding equations for ("!%•) have 
fewer terms than those in Lemmas 12.11 and 1 2 . 2 1 in the fixed metric case. More specifically 
the terms involving the gradient of the Ricci curvature drops out. This is due to a 
cancellation introduced by the Ricci flow. The equation for (wij) will formally stay the 
same. 



Lemma 3.1. Suppose u is a positive solution to \3.1\) such that < u < A. Set 
f = log(u/A) and 



tJ u(l - /) 

with the matrix V = (vij) representing the 2 form ^(T-f) ^ n ^ oca ^ coordinates. Then, 

2f_ 
1-/ 
|V/| 2 , I 



(-d t + A - -r^Vf ■ V)vi, 



j _ j v ij + _ [—2RkijlUki + RuUji + RjiUu] 



18 QING HAN AND QI S. ZHANG 

Proof. This is similar to that of Lemma [2. 11 The only difference is that the commutation 
formula now is 

—dtUijAuij + 2RkijiUki — RikUjk — RjkUik = 0. 
See pl09 of [6J e.g. □ 

Lemma 3.2. Suppose u is a positive solution to hS.l)) such that < u < A. Set 
f = log(u/A) and 

UiUj 

Wij = «2(l_/)2 

with the matrix W = (wij) representing the 2 form u ^^fjy2 in local coordinates. Then, 

v_ 

1-/ 

2|V/| 2 



(-d t +A - -^vf ■ v) Wij 



-Wij + 2(v ki + fw k i)(v k j + fw k j) + RikWkj + RjkWki- 



1-/ 

Proof. This is formally identical to that of Lemma 12.21 The reason is that W = (u>ij) 

du®du — j j„. j-„„ „,uu i-U„ j-;™„ J — ;..„j.;..„ M„™„u, 



represents the 2 form ffit®^ anc : d u commutes with the time derivative. Namely 



, du® du \ dd t u (g) du u® dd t u ( 1 



Hence 



(d t u)iUj Ui(d t u)j t a ( 1 



u 2 (l-f) 2 J u 2 (l-f) 2 u 2 {l-f) 2 1 3 l \u 2 {\- ff 

This is identical to the corresponding term in Lemma [2.21 The computation for all other 
terms are the same also. □ 
In this section we will work with space time cubes that evolve with time. Recall the 
following notation. Let (a?o,io) ^ e a s P ace time point. For R > and T > 0, we write 

QR,T(xo,t ) = {(x,t) | d(x ,x,t) < R, t - T < t < t }, 

We now prove Theorem 11.21 

Proof of Theorem 1 1. 6 A Part (a). Again we set w = tr(W), and 

(3.2) L = -d t + A - r^-jVf • V. 

According to Lemmas 13.11 and 13.21 the following equalities hold 

(3.3) LV = (1 - f)wV + P, 

(3.4) LW = 2(1 - f)wW + 2(V + fW) 2 + Q, 
where P and Q are matrices whose (i,j)-th components are given by 

Pij = 7Z _ 7n [—2Rkijl u M + Ril u jl + Rjl^il] 

(3.5) uyi j ) 

= -ZRkijiVki + Riivji + Rjivu, 
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and 

(3.6) Qij = RikWkj + RjkWki- 
For a constant a G R to be determined, we have 

L(aV + W) = a(l - f)wV + 2(1 - + 2(U + fW) 2 + aP + Q. 

Pick p £ M and a time f where the Ricci flow is defined. Let £ G T p M be a unit 
tangent vector at the point p. Under the metric g(t), we use parallel translation along 
geodesies emanating from p to extend £ to a smooth vector field in the local coordinate 
neighborhood. Then we extend the vector field in time trivially by making it a constant 
vector field in time. We still denote the vector field by £. Since V and W are (2, 0)-tensor 
fields (2- forms), the function 

A = f(aV + W)i = (aV + 

is a well-defined smooth function in a space-time neighborhood of (p,t). Since £ is a 
parallel vector field at time t, it holds, at this time t, 

LX = H + t T (aP + Q)£. 

Here 

H = a(l — f)wfV£ + 2(1 - f)w£ T W£ + 2\(V + /VF)£| 2 . 
By a£ T F£ = A - i T W^, we have 

£T = (1 - f)w(\ - fWS) + 2(1 - f)w£ T W£ + 2\{V + /W)£| 2 

= (i - /)wA + (i - /X T w^ + 2|(y + fw)i\ 2 . 

To simplify the last term further, we fix the space time point (p, t) and assume £ is 
chosen as follows. Under the metric git), we let £ be the time independent vector field 
generated via parallel translation through geodesies emanating from p, by an eigenvector 
of aV + W at (p, t), i.e., at (p, t), 

(aV + W)i = A£. 

Then by (|2,7p in the derivation of (|2,8p . we find that, if A > 0, then at the point (p,t), 
it holds, for a > 4 

Consequently, at time t, it holds 

2A 2 

(3.7) LX > —5- + i T {aP + Q)£. 

Let r be a universal constant to be fixed later. With a = 4, suppose the 2-form 

aV + W - -g(t) 

assumes the eigenvalue /jl = A — f- at (jp\ , ti) , which is the largest for all x G M and 
t G (0, T]. Then A is an eigenvalue of aU + W. Under the metric g(t±), let £ be a 
unit eigenvector of aV + W, which corresponds to A. Under g(ti) again, we use parallel 
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translation along geodesies emanating from p\ to extend £ to a smooth vector field in 
a neighborhood of p\. We still denote it by £. Now we regard £ as a time independent 
vector field defined in a space time neighborhood of (pi,ti). Set 

A = ^(aV + W)i = (aV + W)(£, £)• 



Then, in the space-time neighborhood where £ is defined, we have 

lg(t)(S,S))=L\-L gm ,Z)-l 



We now evaluate at (pi, ii). Since - — !?|?!'^ has its nonnegative maximum at (pi, ti), 



we have, 

(3.9) A(A-^(t)(e,e)) <o, v(A-^(t)(e,o) =o 



> 0. 



and 

r ( a y + ^-ig(t)) 

Here we point out that even though £ is a time independent vector field, its norm changes 
under g(t). This is the reason why we need to normalize its norm in the above inequality. 
Therefore 

d,[(aV + W-fg(t)) ((,()] , r x 2ffeK.fi 
__ + ( aV + W _ _ 9(t) j (tOpj^ £ o. 

Since the computation is at (pi,x±), this implies 

d t (\ - o) + (a - 0)2^, o > o. 



Substituting this together with (|3.2|) and (|3.9p into the left hand side of (|3.8p . we find 
that 

(A - ^)2ffic(£,£) > LX — ^g(t)(Z,0 - -mici^i). 
By (|3.7p . this induces 

2A 2 t 



A2i?ic(£,£) > r-^ - _ - ^(aP + Q)£| at (pi,ii), 



or 



2A 2 r 

(3.10) _ < _ + |^( a p + Q)£|+2A|ffic(£,e)| at(pi,ti). 

To bound A from above, we need to find an upper bound for |£ T \aP + Q)£\ at (pi,*i). 
Let £ = (£i, ...,£ n ). By ()3.5[) and (j3.6|) . we obtain 

|e T (aP + Q)£| < H T Pt\ + \Z T QZ\ 

< |£i£ja (-2R kij iv k [ + RuVji + + |£i£j(-RjfcW fei + R jk w ki )\ 

< [-2R kijl v kl + RuVji + + C|ffic||W|. 
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Writing avki = ctVki + Wkl — Wkl etc in the last line, we deduce 

\a T ( a p + Q)t\ 

(3.11) < [-2R kij i(av k i + wm) + Ru(avji + wji) + Rji(avu + wu)]\ 

+ [-2R m w kl + Ruwji + Rjiwu]\ + C\Ric\\W\. 

At the point (pi,ii), we can choose a coordinate system so that the matrix aV + W is 
diagonal. Let /Ui, • • ■ ,/U n be the eigenvalues of the matrix + — listed in the 
increasing order. We claim that the absolute value of \i\ is bounded from above by fi n 
plus a controlled quantity. Without loss of generality, we assume \x\ < and fj, n = fi > 0. 
Then 

T T 

|^jkijz(a-yfc« +wa:z)| < \Rkiji{av k i + w k i - + \Rkiji~$ki\ 



< | Rkijkjavkk + Wkk - nj) \ + C\Rm\~ 

k=l 



< C\Rm\(ji» + |mi| + j). 

Similarly 

T 

\Rii(ocaVji + Wji)\ < C\Ric\(fj, n + \fii\ + -). 
Combining the last three inequalities, we deduce 

|£ T (aP + Q)£| < C|i?m|(^ n + | M i| + j) + a|i2m||W|. 
In the following, we set K\ = |i?m|i°°. Then 

|e T («P + Q)i\ < CK^n + + \) + C#i|W| 
< Cifi(^ n + | W | + I) + ClfilWl. 

Note 

/il + (n - l)yU n > /ii H h /in 

Hence, 

aA-u r 

Mi < (n - l)/i„ 7- 7T + n-. 

u(l-/) f 

Then, 

\i T {aP + Q)£| < CK X L n - + ^ + C#i|W|. 

Now we need a version of the Li-Yau gradient estimate for the heat equation coupled 
with the Ricci flow. Indeed, by Theorem 2.7 in pQ, we have 

|Vn| 2 u t ci 

2— < — + c 2 K , 



u 2 u t 
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where |-ffic|oo < Kq for some Kq > 0. With ut = An, we get 

Au c\ 

< -f + c 2 K . 

u t 

Since < u < A, we have 

1 

Therefore, at (pi,£i), 

|£ T (aP + < Ci^i (^n + + Cifiif + CK X \W\. 

By the definition of W = (tfy), we have 

iVnl 2 



W < 



u2(l-/)2- 

By [T7j and also [1], we have a curvature independent bound 



\W\<C- f . 



Hence 



\£ T {aP + Q)£| < CK X Un + ^y^J + C^Kq + #i) 

< CK X [\ + i±Ij + C[K y K Q + i^), 

where we used the relation fj, n = ^ = \ — j- < X. Substituting this in (|3.10j) . we arrive 
at 

2A 2 r / l + r\ 

< -2+CK x l\ + ^—\ + C(K 1 K + K 1 ) + 2\K at (pj,^). 

A simple application of the Cauchy inequality yields 



- < 2 Vl + T +ff atfa.ti), 
a t 

where 5 is a nonnegative constant depending only on Kq,K\ and r, with the property 
that B = if K = K x = 0. Then 

T 1 

A - - < (2oVl + T-r)- + aB <aB at (pi,*i), 

by choosing r sufficiently large. In fact, for a = 4, we can take r = 8 + 2\/l7. Recall 
that fj, = A — j at (p±, t±) is the largest eigenvalue of the 2 form aV + W — jg(t) on 
M x (0, T], Therefore, given any nonzero tangent vector rj G T X M, x £ M, it holds 

7/ T (al/ + W)rj- ^g(t)(r],ri) 



g(t)(v,v) 

Thus 

tr] T (aV + WO?? 



<aB inMx(0,T). 



< r + aSt. 



UPPER BOUND FOR HESSIAN MATRICES 



23 



Hence 

t ' \ <- + Bt. 
This proves part (a) of the theorem. 

Part (b). Let ip be a cutoff function supported in the space-time cube Q_r,t(^o>£o) 
such that ^ = 1 in the cube of half the size Qr/2 t/2( x o^o)- We also require that 

„ 12 ) ,v„ £ |, ,a,„c^±1 M,^, EC^. 

Here Kq is a bound on \Ric\ as before. Also the quotients are regarded as when ip = 
somewhere. It is well known that such a cutoff function exists. See [I] e.g. Without loss 
of generality we just take to = T. We also require that -0 is supported in the slightly 
shorter space time cube Qr^t/^{. x Gi *q)- The cut off function can be constructed from 
the distance function. Since we will be differentiating at a fixed point in space time 
eventually, we can use the well know trick by Calabi to get around singular points of the 
distance. 

Now we consider, for some constant f3 £ 1R + to be determined, the 2-form 

^(aV + W)+Pfg(t). 

Let \x be an eigenvalue and £ be a corresponding eigenvector of ip(aV + W) + /3fg(t) at 
some point (p,t), then 

m<*v +w)+ pfgim = rt- 

Here g(t)£ stands for the dual vector of the one form g(t) (•,£,). In a local coordinate, the 
above becomes 

If ip(p, t) ^ 0, then £ is also an eigenvector of aV + W, corresponding to the eigenvalue 
A. Here we just define A by 

H = 0A + Pf. 

As in Part (a), we extend £ to a time independent vector field in a space-time neigh- 
borhood by parallel transport, which is still denoted by £. Now we extend /i and A to 
smooth functions in the same neighborhood by the relation 

l i = fty{aV + W) + Pfg{t)% 

and 

A = £ T (aV + W)£. 
Therefore as functions, fi and A are also related by 

fi = 0A + pf. 

We observe that at the point (p, t), fi and A are eigenvalues of the respective 2-forms. 
However, at different points, this may not be the case. 

Following the computation in deriving ()2.14j) . we know that 

(3.13) 4>L lf x = 0Li(VA + Pf) = i?H + i?i T (aP + Q)£ + ipXLiip + P^L x f. 
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Here L\ is the operator given in (|2.13j) . H is given by (12. 5h and P and Q are given by 
(|3.5|) and (|3.6p respectively. 

Let fl be the parabolic cube given by Q = QR,T(xo,to). Then 



A* 



< 0. 



Here <9 P stands for the parabolic boundary. We will estimate \x from above. 

Since the term H does not involve time derivative, the inequality (|2.6p in the fixed 
metric case still stands. Therefore we have 

?p 2 H > 4t(V>A) 2 + ^ 2 X[w- ArfWn - f^ 2 \ (w- -£ T W£ 
or \ cr / \ a 

At points where /U > 0, we have 

if>\ + Pf>Q, 
and hence ip\ > 0. Then, for a > 4, it holds 

V> 2 # > -^A) 2 . 



This, (I2T51) and (I3T3]) yields 



or 



^iv/i 2 



c 



2|\7^l 2 r- W>' 
IW] + + -1^- + 2y^| V/| • 



V>A. 



As in the previous section, we take 

(3.14) /3 = csup 



with c being sufficiently large. Then we can use the Cauchy inequality to prove: 
(3.15) 



i/>L lfi > \^\? + v 2 e T («p + q)z + l^\vf\ 2 

cr 8 



- C 



\d t i)\ + |A^| + 



C sup 



V> 2 



Let //o be a maximal eigenvalue of ip{aV + W) + /3fg(t) in 0, which associates with a 
unit eigenvector £. Assume /xo is taken at the space-time point (pi,t%). We are interested 
only in the case (jlq > 0. Just like in Part (a), under g(t±), we use parallel translation 
along geodesies emanating from p\ to extend £ to a smooth vector field in a neighborhood 
of pi. We still denote it by £. Now we regard £ as a time independent vector field defined 
in a space time neighborhood of (pi,ii). Set 

ll = t T [^(aV + W) + Pfg(t)]i. 



Since 



(3.16) 



has its nonnegative maximum at (pi,ii), we have, at this point, 
An = &(n/g(t)(£, 0) < 0, Vm = V(/z/s(t)(& 0) = 0, 
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and 



Therefore 



dt 



> 0. 



g(m,0 r \g(m,0\ 2 

Since the computation is at (px,xx), this implies 

d tf jL + fJ,2Ric(£, > 

Recall from (j27T3j) that 



> 0. 



L a = -d t + A 



2/ ; V/-V-^V. 



1-/ V 
Hence we can plug the above inequality and A3. 16H in (|3.15p to deduce 



1 



C 



|^| + |A^| + 



C sup 



This implies, at 



— o(<M) < r ie J (aP + Q)£| + C - + + rl>n2K . 



mi 



(3.17) a ,„„, VT . ^ 

Now we control the right hand side. From (|3.1ip in Part (a), we have 
^\f{aP + Q)i\ 

< [-2Rkijiip{av k i + to^) + Rui>{avji + w^) + Rjiip(av u + 
+ ^ [-2R kij iw M + R u wji + flj-jwall + Cty|12ic||W| 

< [-2Rkijii>(av k i + w H ) + Ruip(avji + ii;^) + Rjiip(avu + ?%)]| 
+ (7(^0 + ^)^1 

< [ - 2R kij i{^(av kl + u> fc/ ) - /3/5 fc /} + Rii{^(avji + w jt - f3f8 H } 

+ Rji{ip(avu + u> i; - /3/5fcj}] | 

+ C(lfo + K x )Pi>\f\ + C(K + ETi)^|iy|. 

Let Hi, ...,/u n be the eigenvalues of the 2-form ip(aV + W) + /3fg(t) at the space-time 
point (pi,t\), which are listed in increasing order. We assume without loss of generality 
that [i\ < 0. Then the above inequality implies 

iP\t T (aP + Q)£| < C{K + K x ){ii n + \^\) + C(K + + C(ET + K{)^\W\. 

Observe that 

/ii + (n - l)/i n > jtii H h 



tr 



an. 



+ 



u 2 (l-/) 2 



'7 



u{l ~ f) 
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Hence, 



|MX I < (n - l)A* n - 7T + 

«(1 - /) 



Therefore, 



(3.18) Mf{aP + Q)£| < CiTi ( Mn - V'^^tt ) + CK^\W\ + CK x fi\f\. 



u(l ~ f) 



By Theorem 2.7 in [T] again, we have 



, |Vu| 2 ut\ ( 1 1 . 



u- 



l + t R 2 



where Kq = \Ric\oo- We mention that this inequality is not exactly the one stated in 
their Theorem 2.7, due to the appearance of the cutoff function i/j in the front and the 
appearance of instead of ~. However, this is actually what was first proved at the 
first line in p3532 there. This is also the case for the Li-Yau inequality in the fixed 
metric case. Since, by construction, ip is supported in the slightly shorter space time 
cube Qr^t/a ( x o > to) an d to = T, we have 

Using ut = Au, we get 

, ( \ n 

-ip — < c 3 - + + c 2 K . 



u ~ \T R 2 

Since < u < A, we have 

1 



< 1. 



1-/ 

Substituting these in f|3. 18|) . we find that, at (pi,ii), it holds 

il>\f{aP + Q)£| < CK X (\ + c 3 (~ + ±\ + c 2 K^j + CK^\W\ + CK x p\f\. 

Here we also used the inequality \x n = \i = ip\+0f < ip\. By the definition of W = (wij), 
we have 

According to Theorem 2.2 in PQ, we obtain 

Again the cutoff function ij; and term ^-j- were not in the original statement. But this 
was proven in that paper on the way to prove Theorem 2.2 there. Also the last inequality 
follows from the choice of the support for ip. Now we know that 

^ 2 \e(aP + Q)£| < CK^X + CK^ + C{K X + 1)K + CK x ^\f\. 
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Substituting this in (|3.17j) . we conclude, at (pi,t\) 



1 o „„ .. „„ / 1 1 \ „ f\ 1 X 2 



< CK^X + CR\K + CT^ M + _ + _ j + + _ 

+ CK x ^\f\ + 2K^ix. 
As n = ipX + /3/ < ip\ and K < CK U this yields 



^x<c(^ + ^ + b\+c^k^\T\, 



where B is a nonnegative constant depending only on Kq,Ki such that B = when 
Kq = K\ = 0. This shows 

n = <M + pf < c (± + ^ + b\ + cVifi Wl + Pf- 



Since / < 0, we know that Cy/Ki/3tp\f\ + /?/ < CKi, which implies, at (pi,ii), that 

ti = r/>\ + Pf<c(± + -^ + B 

Recall /i at (pi,ti) is the maximum of the eigenvalues of the 2-form ip(aV + W) + f3f 
in the cube Qr,t{xq, to). Hence for any unit tangent vector £ at a space-time point 
(x,t) G QR,T(%a,to), it holds 

« T K + W0£ + /?/ < c ( - + 4 + 5 1 > 



T i? 2 
or 

^e T («^ + < c ( - + i + + /3|/|. 

With the choice of /3 = csup = ^ in (|3,14p with c sufficiently large, we then have 

^ t v^<c(^ + ^+b\ (1-/), 

and hence 

^ £C (? + ^ + B ) (1 - /)2 - 

This implies the desired estimate. 

□ 
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